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ABSTRACT 


As a part of the “Dust, Ice, and Gas In Time” (DIGIT) key program on 
Herschel, we observed GSS30-IRS1, a Class I protostar located in Ophiuchus 
(d = 120 pc), with Herschel /Photodetector Array Camera and Spectrometer 
(PACS). More than 70 lines were detected within a wavelength range from 50 
/mi to 200 /mi, including CO, H 2 0, OH, and two atomic [O I] lines at 63 and 
145 /mi. The [C II] line, known as a tracer of externally heated gas by the 
interstellar radiation field, is also detected at 158 /mi. All lines, except [O I] 
and [C II], are detected only at the central spaxel of 9".4 x 9".4. The [O I] 
emissions are extended along a NE-SW orientation, and the [C II] line is detected 
over all spaxels, indicative of external PDR. The total [C II] intensity around 
GSS30 reveals that the far-ultraviolet radiation field is in the range of 3 to 20 
Co, where Go is in units of the Habing Field, 1.6 x 10~ 3 erg cm -2 s -1 . This 
enhanced external radiation field heats the envelope of GSS30-IRS1, causing the 
continuum emission to be extended, unlike the molecular emission. The best- 
fit continuum model of GSS30-IRS1 with the physical structure including flared 
disk, envelope, and outflow shows that the internal luminosity is 10 L 0 , and 
the region is externally heated by a radiation field enhanced by a factor of 130 
compared to the standard local interstellar radiation field. 


Subject headings: 
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INTRODUCTION 


The shocked gas induced by interaction between the protostellar jet and the surrounding 
interstellar medium cools through copious molecular emission lines and forbidden atomic 


lines, s uch as CO, H 9 Q. OH , and [O I] lines, in the far infrared (FIR) range (INisini et al 


2002|). iGiannini et al.l (120011) suggested that the FIR lines are good tracers of the energy 
budget related to outflows based on the observations of young stellar objects (YSOs) with 
the Long Wavelength Spectrometer (LWS) on board the Infrared Satellite Observatory 
(ISO). According to their results, the powerful outflows of Class 0 sources are the results 
of larger envelope mass infall rates as compared with more evolved sources, indicating that 
outflow can be an evolutionary parameter of YSOs. Therefore, the observations in the FIR 
band are important to understand the nature of the early stage YSOs. 

The [O I] and [C II] lines at FIR wavelengths have been used to measure the strength of 
the interstellar radiation field (ISRF) in order to predict the presence of photo-dissociation 

and the heating mechanism associated with YSOs 


regions (P 

DRs) as well as underst 

(Yui et al. 

1993; 

Nisini et al. 

2002) 


2002). Although the enhanced interstellar radiation fields 


result in the high excitation conditions, the effects of the external radiation fields on the 
low mass YSOs are poorly understood. 

Observations of YSOs using ISO /LWS stimulated studies of origins of the FIR emission 


as w ell as line cooling budgets in atomic and molecular line emission (lOiannini et al 


2001 


Nisini et al. 


2002). However, the ISO beam was too big to discriminate the different gas 


conditions. Recent 


2013 : 


Karska et al. 


, lar ge numbers of embedded YSOs have b een observed (Green et al. 


2013) with the Herschel Space Observatory (IPilbratt et al. 


2010) with 


much better spatial and spectral resolutions and much better sensitivity. Many weak 
molecular and atomic lines have been det ected with the Photodetector Array Camera 


and Spectrometer (PACS; 


Poglitsch et al. 


2010) on board Herschel. Here we present the 
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detailed analysis of line and continuum fluxes of GSS30-IRS1, which was observed by the 
Herschelj PACS instrument as part of the Herschel Key Program, DIGIT, “Dust, Ice, and 
Gas in Time” (PI: N.Evans). 


The Ophiuchus 


(iLoi nard e t ah 


(hereafter Oph) molecular cloud, at a distance of 120 ± 4 pc 


2008 1. is a good region for studying the impact of irradiation fields because 


the interstellar radiation fiel d impinging on th e Oph cloud is st r onger than the average 


for our solar neighborhood (jLiseau et al 


1999). 


Lindberg et ah 


(2014) showed that the 


external radiation field could affect the physical and chemical properties around YSOs 
based on Herschel observations. 


One of the well studied targ ets, GSS30, is a bipolar reflection nebula in the Oph 


molecular cloud. 


Weintraub et ah 


(119931 1 reported that there are three infrared sources: 
GSS30-IRS1, -IRS2, and -IRS3 (Figure Hh£ According to the geometry of the infrared 


emission, GSS30-IRS 


in the near-infrared ( Oastelaz et ah 


is the illuminating source in the nebula and the strongest sou rce 


1985; 


Weintraub et ah 


1993 


Zhang et ah 


1997) 


(Figure dh>). GSS30 


-IRS1 has been classified as a C l ass I YSO by its sp ectral energy 


distribution (SED) (IGrasdalen et a 


L q and Tboi of 138 K (1 Green et al 


1973 


Elias 


1978 


Wilking et ah 


1989 ) with L bol of 14.5 


2013). IRS2 has no detected emission at millimeter 


wavelength, indicating that it is_a more evo 


IRS3 was detected at 2.7 mm (Zhang et al 


19971 and 6 cm (ILeous et al 


ved (T-Tauri) star f An dre fc Montmerle 


known as a Class I YSO with L b oi of 0.13 L Q (IBontemps et al 
IRS1. 


1994) 


1991)- It has been 


20011) . much fainter than 


The origin of the associated molecular outflow is unclear. 


Tarnura et ah 


( 19901 ) detected 


high velocity millimeter CO 1-0 line emission in this region of roughly 5' in extent. However, 


1 The IRAC image was downloaded from the Spitzer data archive. 
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both the blue- and red-shifted lobes lie to the south of GSS30-IRS1, deviating from the 
bipolar morphology seen in the reflection nebula. They interpreted this asymmetry of the 
CO outflow lobes in terms of the interaction between one side of the lobe and a dense NH 3 
cloud core which is located to the northeast of GSS30-IRS1. 
that high velocity gas components reported by 


Andre et al. 


Tamura et al 


the VLA1623 jet, which is located southeast of GSS30. Recently, 


(11990!) suggested 


1990 ) are associ ated with 


Jprgensen et al. 


(1200911 


detected wing emission clearly associated with the GSS30-IRS1 in the HCO + lines observed 
with the Submillimeter Array (SMA). 


In order to study the effect of external radiation fields on GSS30-IRS1, we analyze the 
observed PACS atomic line spectra and model the continuum emission using the continuum 
radiative transfer code (RADMC-3D). For molecular lines associated directly with the 
protostar, the rotational diagrams of molecular transitions and a non-LTE large velocity 
gradient (LVG) code (RADEX) have been utilized. The outline of this paper is as follows: 
in § El we describe the Herschel observations and data reduction. The observational results 
are presented in § El The detailed analyses, including line and the continuum emission, 
are described in § |4] and § El We also discuss the effects of external heating in § El and 
summarize our results in § El respectively. 


2. OBSERVATIONS AND DATA REDUCTION 

GSS30-IRS1 was observed on 10 March 2011 (A ~ 50-75 and 100-145 /mi; AOR: 
1342215678, A ~70-100 and 140-210 /mi; AOR: 1342215679) by Herschel /PACS in a 
pointed range-scan spectroscopy mode. The PACS instrument consists of a 5 x 5 array 
whose unit spatial pixel (hereafter spaxel) covers 9".4 x 9".4 with a spatial resolution of ~ 
9" at the shorter wavelengths and ~ 18" at longer wavelengths. Its spectral coverage ranges 
from 50 /mi to 210 /mi with a spectral resolution of X/S\ ~ 1000-3000, and its coverage is 
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divided into four segments, covering A ~ 50-75, 70-105, 100-145, and 140-210 /mi. The 
edges of the segments (A < 55 /mi and 95 < A < 103 /mi) suffer from noise, and the spectral 
fluxes appear much higher than their real fluxes longwards of 190 /mi. For those reasons, 
we limited the wavelength ranges to 55 < A < 95 /mi and 103 < A < 190 /mi in this study 
to calculate the line fluxes accurately. The sky background and telescope emission were 
subtracted using two nod positions located 6' from the source in opposite directions. 


The basic data reduction of the DIGIT obse rvat i ons w as performed with Herschel 


Interactive Processing Environment (HIPE) v.13 flOtt 


2010). The complete description for 


the data reduction can be found in Green et al. (in prep.). This reduction provided the 
best signal-to-noise ratio (S/N) and better absolute flux calibration. In order to measure 
the total line fluxes, we used the two different methods due to the spatial distributions of 
molecular and atomic lines. For all molecular lines, the spectra extracted from the central 
spaxel were corrected for the point spread function (PSF) of PACS. However, we did not 
apply this method to atomic lines because the [O I] and [C II] lines are detected over all 25 
spaxels. In addition, the PACS fluxes for the atomic lines are contaminated by the emission 
in the sky positions (Sect. 13.4p . Thus, we measured the fluxes from the nod B observation 
(with off position 6' to the south from GSS30-IRS1), in which [O I] and [C II] appear in 
emission. While the [C II] lines are analyzed in all 25 PACS spaxels, the [O I] fluxes are 
measured over only 15 spaxels (red crosses in Figured]) selected to avoid the effects of the 
outflow by VLA1623 in the southwest spaxels (see Sect. 13.31 and Sect. 13.4[) . 


All lines were fittec 
(SMART; Higdon et al. 


by the Spectroscopic Modeling Analysis and Reduction Tool 


200 4!). and the Leiden Atomic and Molecular Database (LAMDA) 


was used as the molecular database flSchoier et al. 


2005). 
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3. RESULTS 


3.1. Spatial Distribution of Continuum and Line Emission 


One interesting observational result is that the continuum emission extends to the 


northeast. The same feature has 
the Very Large Telescope ( VLT ) 


leen shown as sc attered light in an H -band image from 


Chen et al 


20071 ) and a 1.6 /mi Hubble Space Telescope 


(HST) image from the Near Infrared Camera a nd Multi-Object Spectrom eter (NICMOS) 


(Allen ct ah 


200 2). The SCUBA 850 pm map ( Pi Francesco et al. 2008 ) also shows the 


same extended structure (Figure [T}i). In Figure [21 the continuum emission is significantly 
more extended than the line emission at longer wavelengths except for [O I], although a few 
lines at longer wavelengths, A > 120 /mi, appear extended, which is mostly because of the 
large PACS PSF. Thus, we explore the cause of the extended dust continuum emission by 
modeling the SEP of GSS30-IRS1 in Sect. El 


3.2. Continuum Deconvolution 


As seen in Figure Ek,, there are three YSOs within the PACS field of view (FOV). 
Thus, it is difficult to measure the total flux of IRS1 from the PACS total fluxes because 
of the contamination by those multiple sources. In addition, the PACS image at 70 /mi 
(Figure [T}d) shows that the continuum is clearly extended. In order to test whether this 
extended feature is due to the contribution from three YSOs or the external heating, we 


decompose the fluxes of IRS1, IRS2, and IRS3 with the assumption that a 
point-like. The key concept of deconvolution method is the same as used in 


1 sources are 


Lindberg et ah 


( 20141 ). except that we adopted the PSF from the PACS spectrometer beams version 


2 See https://nhscsci.ipac.caltech.edu/sc/index.php/Pacs/PACSSpectrometerBeams 






















Figure [3] shows the deconvolved spectra of three sources in the PACS range. All spectra 
are deconvolved within the PACS FOV. Tims, we note that deconvolved flux densities of 
IRS2 and IRS3 are underestimated because they are located near the edge of the PACS 
FOV. However, the deconvolved flux densities of IRS1 are similar to the fluxes measured 
from the central spaxel and then corrected for PACS/PSF, within the calibration error (10 
%). It indicates that IRS1 is located at the central position of PACS and the deconvolved 
flux densities are the total flux densities of IRS1. The continuum emission from IRS3 is 
comparable with that from IRS1 beyond 100 pm (Figure [3]). 

Figure H] shows the residual emission after subtracting emission by individual infrared 
sources from the total fluxes over 5x5 spaxels as described in the equation below: 

Fv, residual F V) 5 x 5 (-^1/, IRS1, deconv. T F u , IRS2, deconv. T F v ^ IRS3, deconv.)) (1) 

F Uy 5 x 5 is the flux measured over the PACS 25 spaxels and F v narn p deconv. are the 
deconvolved fluxes of each source in the PACS FOV. The residual spectrum is fitted by the 
black-body temperature of ~ 40 K. 


3.3. The Spectrum of GSS30-IRS1 


In Figure [51 we present the SED of GSS30-IRS1 up to 1300 pm obtained from the 
photometry of 2MASS and Smtzer -IRAC and the spectroscopy of Herschel /PACS and 


Spitzer -IRS (iLahuis et ah 


2006) from the Spitzer Legacy Program “From Molecular Cores 


to Planet Forming Disks” (c2d) (Evans et al. 


Time Observations (GTO) program ( Furlan et al. 


20031 and the “IRS_Disk” Guaranteed 


2006). The sub-mm data are from 


the Submillimeter High Angular Resolutio n Camera II (SHARC-I 


Common-User Bolometer Array (SCUBA) ( 


van Kempen et al. 


Radioastronomie Millimctrique (IRAM) ( Motte et al.l 1998). The deconvolved SED of IRS1 


), Submillimetre 


2009), and the Institut de 
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(Figure El red line) exhibits a bolometric luminosity of 8.5 L & and a bolometric temperature 
of 192.8 K. 


The rich PACS spectra show molecular lines including CO, H 2 O, and OH, and two 
atomic [O I] lines in emission. However, as seen in Figure El the [C 11] 158 /ini line has a 
deep absorption feature. All detected lines are listed in Table Q] with blended line emissions 
noted. We excluded blended lines in our analyses. Detected lines are presented in Figure [6] 
to Figure [9] Although weak CO lines at long wavelengths are detected toward spaxels 
south of the central source, those lines probably trace the outflow related to a nearby 


Class 0 source, VLA1623, located southeast (5" offset in rig’ 


it ascension and V offset 


in declination from GSS 30-IRS1) of GSS30-IRS1 ( Dent et al. 


Pontoppidan et al. 


1995; 


Yu & Chernin 


1997; 


2002). Blue and red outflow win gs, re late d to VLA1623 and other HH 


objects, are clearly seen south of GSS30-IRS1 (IBierkeli et al. 


2012). In addition, the 12 CO 


J = 3-2 observati ons s how the outflow emission dominated by VLA1623 toward the south 


of GSS30-IRS1 (van der Marel et al. 


2013). To avoid contamination, we present only the 


CO, H 2 0, and OH lines extracted from the central spaxel (Figure El to Figure [9]). 


In contrast, two detected atomic lines ([O I] and [C II]) are distributed over a region 
larger than the central spaxel and they are often affected by strong emission in the nod 
positions used for sky subtraction. We investigate the spatial distributions of these lines 
using two different sky subtractions in the next section. 


3.4. Spatial Distribution of Atomic Lines 

Figure fTITl shows the entire PACS spectral line map of the [O I] 63 /im line, which is the 
average of the two nods subtraction. The [O I] line emission is elongated along the NE-SW 
direction, and the morphology is similar to a part of the outflow cone detected by the SMA 
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HCO + 3-2 and HCN 3-2 observations (Jprgensen et al. 


2009). This indicates that the [O 


I] traces the outflow structure (see next section). Although [O I] emission extends into the 
region where IRS3 is located, we conclude that the [O I] emission is related to IRS1, due 
to the strong peak at the IRS1 position in our line map. In ad dition, the outflow tracer 


HCO + and HCN lines are clearly detected toward IRS1 ( Jprgensen et al. 2009). There is a 


separate structure to the SW of GSS30-IRS1 in the line map of the [O I] (Figure fTUl). As 
mentioned above, the structure seems related to the outflow by VLA1623. 


Figure fill to Figure ITHI show the atomic line maps where the spectrum of each nod 
position was subtracted separately. The [O I] line at 63 /iin is clearly seen as an emission 
line in both subtractions and its morphology is similar to that found in the final data 
reduction (Figure [Toll . suggesting that it is less affected by emission in reference positions. 
In maps of the [O I] 145 fin l and [C II] 158 fin l (Figure fl2l and Figure fT3l) . however, the left 
panel subtractions show an absorption feature because of the emission in the nod A position 
itself. However, the other subtractions in the right panel, using nod B, show emission lines. 
Therefore, the relative distribution of the [O I] 145 fim and [C II] 158 /im emission using 
nod B could be on-source emission although the actual intensity would be stronger if [O 
I] 145 /jrn and [C II] 158 /mi emission exists in the nod B position. The [C II] intensity 
contour map (Figure fITll uses the nod B subtraction in Figure fl3b. The [C II] emission 
does not peak in the central spaxel and it has no spatial correlation with [O I] emission 
(Figure [2]o and e), indicating that they are not excited by the same mechanism as [O I], 

Although there is a possibility that the atomic lines are contaminated by the nod 
positions, in this analysis we use the atomic lines of the nod B subtraction assuming the 
result is representative of the source. 
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LINE ANALYSIS 


Recently, a number of studies have presented the results of rotational diagrams and 


non-LTE an alysis of molecu lar 


Green et al 


Green et ah 


2013 


Karska et al. 


ine emission detected with Herschel /PACS flNeufe 


2013 


Lee et al. 


2013; 


Manoj et al. 


2013 


Lee et al. 


20 


2014]). 


(120131) has described the simple rotational diagram of all DIGIT sources with 


molecular lines including GSS30-IRS1; GSS30-IRS1 is generally consistent with other 
embedded sources (Green et al. in prep.). 

As clearly seen in Sect. [3j there is a different trend between atomic and molecular lines; 
molecular lines are centered on GSS30-IRS1 and are more compact compared to atomic 


lines. In the Taurus sources, [C II] emission is to o weak to see t 


most sources, except for L1551-IRS1 and TMC1 (Lee et ah 


le spatial distribution in 


2014). The [C II] emission in 


those two sources is well-correlated with [O I] emission, related to the outflow direction. 
However, the spatial distributions of [O I] and [C II] lines in GSS30-IRS1 have no correlation 
with each other. We provide the results of the rotational diagram and non-LTE analysis 
of molecular emission without much detailed description because the results are not very 
different from previous work toward other embedded YSOs. In this paper, we will focus 
more on the analysis of atomic lines. 


4.1. Molecular Lines 


Rotational Diagrams 


In the rotational diagram, CO ladde rs are fittec 
breakpoint at E up = 1800 K as done in 


Green et ah 


yy tw o rot ational temperat ures with a 


( 20131 ) and 


Lee et al. 


(120141 . We refer 


to those papers for the fundamental concept of rotational diagrams. For the “warm” CO 
component, the rotational temperature is T rot ~ 359 ± 4 K, and the “hot” component is 
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fitted by T rot ~ 783 ± 32 K. However, the apparent rotational temperatures of H 2 0 (T rot 
~ 197 ± 2 K for ortho-H 2 0, and T mt ~ 171 ± 5 K for para-H 2 0) and OH (T rot ~ 193 ± 7 
K for 2 n 3/2 ladders, and T rot ~ 128 ± 5 K for 2 Tii /2 ladders) are much lower than that of 
CO. H 2 0 and OH are usually sub-thermally excited due to their high critical densities. The 
derived rotational temperatures and total number of each molecules agree well with other 
DIGIT embedded sources. All results of rotational diagrams are found in Table [21 


Non-LTE Analysis 


For a more realistic : analysis o: 


probability code (van der Tak et al 


molecular excitation, we utilized the RADEX escape 


2007) to understand the physical conditions surrounding 


GSS30-IRS1 assuming a plane-parallel molecular cloud. We adopted the sam e ranges of 


physical parameters (T kin , n(H 2 ), and N (molecule)/A v) as used in 


Lee et al. 


(2014) to find 


the best-fit models. 


Lee et al 


(120131 ) and 


h e det ails of the modeling procedure with RADEX can be found in 


Lee et al. 


(2014). According to the non-LTE LVG models, the CO 


emission lines favor the highest temperature (7km = 5000 K), and a power-law temperature 
distribution can fit the CO emission well, indicating the sub-thermal solution of CO. 


Our results for CO show a similar trend to those found in TMC1 and TMR1 (L ee et ah 


2 0141 ). For H 2 0, one component model explains the observed H 2 0 emission better than the 
power-law model since the power-law model is meaningful only when the lines are optically 


thin. The best-fit model temperature (500 K) in the one component moc 


lower than that of CO or that of H 2 0 for all Taurus sources (Lee et al. 


el for H 9 O is much 


Lee et ah 


2 0141 ). For OH, 


(120131 ) reported that the IR-pumping plays a role only for the highest energy 


level (E up = 875) transition. However, our best-fit model suggests the model without 
IR-pumping fits better the highest energy transition than the model with IR-pumping. The 
model parameters of our best-fit model of each molecule are listed in Table [3j 
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4.2. Atomic Lines 


The Origin of the [C II] emission 


The Oph molecular cloud is a good region to study the [C II] emission, which is known 
as a tracer of the PDRs excited by stellar or interstellar radiation field: (1) the cloud 
contains three B-type stars, (2) is close to us, and (3) is at a high enough galactic latitude 


to avoid the emission from the galactic plane flYui et al. 


1993 ). Our PACS o bserva tions also 


show extended [C II] emission around GSS30-IRS1 (Figure mj). According to 


Giannini et al. 


(120011 1. the detection of [C II] lines originating from embedded low luminosity YSOs is very 
unlikely due to their low luminosity and high density of envelope. Therefore, the detection 
of the [C II] line without any correlation of its spatial distribution with any other molecular 
emission indicates an external PDR activity around GSS30-IRS1. 


Yui et al. 


(1199311 and 


Liseau et al. 


(11999 1 reported three possible candidates of external 
UV sources: HD 147889 (B2 V), S 1 (B3 V), and SR 3 (B7 V) in the Oph cloud. Their 


visual extinctions are A v = 4.6, 12.5, and 6.2 mag, respectively (jElias 


1978). Both studies 


assume that HD 147889, which is the hottest among the three stars, is the dominant 


external U V sou rce, thoug 


h S 


model by 


Kaufman et al 


is closer in projection to the Oph cloud. Based on the 


(119991 ). we interpret the strength of the far-ultraviolet (FUV) 


radiation field, Go (in units of the Habing Field, 1.6 x 10 3 erg cm 2 s 4 ) combined with 
density of H nuclei (n) and [C II] 158 /ini intensity. 


As we measured the [C II] line fluxes from the spaxels where [C II] line is detected in 
emission (Figure U3b). we can calculate the total intensity of [C II] lines using emitting area 
with an assumption that the sky emission at the nod B position is insignificant. 

The intensity of [C II] calculated over the PACS FOV is 1.1 x 10~ 5 erg cm -2 s -1 sr” 1 . 
If we adopt the density of 10 4 < n(H) < 10 5 at the radius greater than 2000 AU, which 
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is calculated from our dust continuum modeling 


range of 3 to 20 based on the PDR model from 


(see Sect. [5]) 2 th e stre ngth of Go is in the 


Kaufman et ali (110001 ). Liseau’s empirical 


model, based on their observed angular distribution of the [C II] line fluxes, also reveals 
that the value of Go in the Oph cloud is in the range of 10 to 100. The derived Go value 
gives us the lower-limit because there might be the [C II] emission in nod B position while 
we assume no emission at the nod B position. 


The Origin of the [OI] emission 


Although not as significantly as the [C II] emission, the [O I] lines at 63 and 145 /irn 
are also affected by the sky emission as seen in Figure fTTk and Figure fl2h. In particular, 
no emission line at 145 /irn is seen at the central position of Figure fl2h . indicating that 
the sky emission from the nod A position is comparable to the source flux. In order to 
calculate the actual intensity of [O I], we applied the same method to [O I] lines as done for 
[C II] except that the [O I] lines are measured from the PACS 15 spaxels (red crosses in 
Figured])- The [O I] intensity is 2.2 x 10 -4 and 1.2 x 10 -5 erg cm -2 s“ 


1 2 ° 1 sr 1 for 63 and 


145 /irn, respectively. The derived [O I] 
the strength of Go can reach up to 10 6 flKaufman et al¬ 


ine ratio between 14 5 /rm and 63 /iin suggests that 


19991) . which is unrealistically high 


compared to other studies and our [C II] analysis. 


The [O I] emission seems associated with the HCO + flow ( Jorgensen et al . 2009), which 


is extended beyond the PACS FOV. It is also peaked at the central position both at 63 and 
145 /irn (Figure [2 Jd and e). This line morphology suggests that the [O I] emission is more 
directly related to the source itself, unlike the [C II] line. Therefore, a fast dissociative 
shock, which is produced by strongly collimated jets, or an irradiated C-shock along the 
outflow cavity wall, which is exposed by the UV photons from the central source, may cause 
most of the [O I] emission in GSS30-IRS1. 
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The interactions between the outflow/jet and the dense envelope produce highly 
energetic shocks which propagate into molecular gas, consisting of H 2 , CO, and H 2 0. The 
abundances of those species strongly depend on the shock condition because the shocks 
could dissociate the molecules. Therefore relative intensities of those gas emission lines 
can characterize the type and speed of associated shocks. Thus, we com pare the observed 


line fluxes with the shock models by Fl o wer fe Pi neau des Foret s! (l2010h . who presented 
models for C- and J-type shocks, where the shock speed, v Sl ranges from 10 km s" 1 to 
40 km s'" 1 , and pre-shock densities, nn, from 2 x 10 4 cm" 3 to 2 x 10 5 cm" 3 . The shock 
speed of 40 km s" 1 is only for the C-shock. According to the [O 1] 63/145 /irn flux ratio in 
GSS30-IRS1, the [O 1] flux ratio is well fitted by a high density C-shock model (Figure fl5l). 


We al so co mpared th e actu a 


intensity of [O I], which is calculated above, as done in 


Flower fe Pineau des Foretsl (120101) . Both [O I] lines are matched by a high density J-shock 


model with = 17 - 27 km s" 1 fFigure [T6lh The [O 1] 63 //m line is excited more easily than 
the [O I] 145 /irn line due to its lower excitation energy (228 K). Therefore, if the J-shocks 
actually play an important role in GSS30-IRS1 as expected by the actual intensities of the 
[O I] lines, the low [O 1] flux ratio, which is better explained by C-shocks, suggests that the 
[O I] 63 /in i line is optically thick. 


4.3. Line Luminosities 


In Table 01 we present the FIR continuum and total line luminosities in the wavelength 
range covered by PACS as well as the FIR line lumin osities o f each species. The definition 


Nisinietah 


( 2002 ): Lfir — Fco + 


of total FIR line luminosity, Lfir, was adopted from 
Lh 2 o + Lqh + Ly oi]. In the PACS wavelength range, fractional contribution of the total 
line luminosity to the total FIR continuum luminosity is ~ 1.9 x 10" 3 , indicating that the 
cooling by the dust is dominant in FIR radiation. 
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According to I Oi ami ini et al.l (120011) and iNisini et al.l (120021) . the ratio of the total 
molecular cooling to the bolometric luminosity (L mo i/L bo i) in the FIR wavelengths is 
significantly higher (h mol /L bol > 1 x 10~ 3 ) for Class 0 sources than for Class I sources 
(hmoi/hboi < 10~ 3 ), suggesting that the outflow power is greater in Class 0 than Class I. In 
GSS30-IRS1, the ratio of the FIR molecular cooling in the PACS range to the bolometric 
luminosity is ~ 8 x 10 -4 , which is consistent with the Class I range. 

The percentage contributions to total line cooling are 21, 36, 11, and 32 % for CO, 
H 2 0, OH, and [O I], respectively, in GSS30-IRS1. Thus the cooling through molecular 
emission is dominant (68 % to Tfir), but [O I] is also a significant coolant (32 %). By 


contrast, analysis of the Class I protostar L1551-IRS5 indicates [O 


70 % to the total line cooling, with almost no contribution of H 2 0 (Lee et ah 2 0141 ). The 


contributes around 


36 % contributed cooling from H 2 0, similar to [O I], in GSS30 is indicative of the lower 
dissociation rate of H 2 0. 

If GSS30-IRS1 is associated with a PDR excited by the UV photons from a nearby 
luminous B type star or accretion to central pro tostar, t he_abundances of [O I] and OH 


would be enhanced compared to that of H 2 0 (iHollenbach fe Tielens 


1997 ) due to the 


dissociation of molecules. We t 


ius calculated 


the ratio s with Taurus sour ces (ILee et al. 


he total line luminosity ratios and compared 


sources (jKar ska et al. 


2014), L1448-MM ( Lee et al. 2013), and all WISH 


2013). We find that GSS30-IRS1 has a small OH/H 2 0 luminosity 


ratio compared to other embedded sources (Figure fTTl) . This does not necessarily imply that 
H 2 0 is not being dissociated, since the OH, the transient species in the dissociation of H 2 0, 


survi ves for a very short time before its dissociation to H and O (Ivan Dishoeck fe Dalgarno 


1983). The [O I]/H 2 0 ratio of GSS30-IRS 


influence of high energy photons (Lee et al 


is c 


2014). 


ose to unity, indicative of the possible 
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DUST RADIATIVE TRANSFER MODELS 


Internal Heating 


Kristensen et al.l (120121) found the envelope physical structure of GSS30-IRS1 by fitt ing 


its SED with the 1-D continuum radiative transfer code DUSTY (Ivezic & Elitzur 


1997). 


They considered a density profile, n oc r~ p , where n is the density, r is the radius, and p 
is the power-law index. The fitted SED data covered from 24 /irn (Spitzer- MIPS) to 850 
/irn (SCUBA) including the Herschelj PACS observations. However, the NIR and MIR 
datapoints cannot be fit by this 1-D envelope, instead requiring a disk structure, which 


radiates predominantly at shorter wavelengths (Kristensen et al 


2012). In addition to the 


disk, embedded YSOs have strong outflows to develop bipolar cavities. Thus, we construct 
a 2-D SED model using the deconvolved spectrum of IRS1 assuming that IRS1 is the 
dominant internal heating source in the vicinity of GSS30. 

To explore the disk and envelope structures of GSS30-IRS1, we modeled the SED 


using the Monte Carlo radiative transfer package RADMC-3D v3.1 (Dullemond fe Dominik 


2004 ). The main inputs are the density structure of the dusty circumstcllar material, the 
dust opacity tables, and the stellar parameters. The model assumes that the s catte ring 
is isotropic. We adopted the same dust grain models used by 


mixture of silicates and graphite grains ( jLaor & Drai ne 


Robitaille et ah 


(2006), a 


1993). We used different grain 


models for different regions in the di sk and t he envel ope a s follows: “Disk midplane” dust 


grain model described in Table 3 of 


Whitney et al. 


( 20031 ) for the dense region (nn 2 > 


0 10 cm 3 ) and 


(Kim et al. 


lie “KMH” dust grain model for the rest of the circumstellar geometry 


1994)- The grain model for the disk midplane has a size distribution which 


decays exponentially for grain sizes between 50 /irn and 1 mm. The grains in the_ envelope 


are s lightly 


(2003) and 


larger than ISM 


Robitaille et ah 


grains. For the detailed grain properties, see 


( 12006) . 


Whi tney et a l. 
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To determine initial parameters, we used the online SED fitting tooj^l by 
( 20061 ). This tool offers a grid of 200,000 SEDs including stars, disks, and envelopes with 
a wide range of physical parameters. The SED is sensitive to the total internal luminosity 
and the inner density of the envelope, especially in FIR range. The total internal luminosity 
is the sum of the stellar luminosity and accretion luminosity, where the stellar luminosity 
is given by L* oc R 2 * T 2 * and accretion luminosity, L acci is a function of stellar parameters 
and disk accretion rate Md; s k- Since there is high degeneracy between accretion and 
stellar luminosity, we simplify the fitting by adjusting the total luminosity only via stellar 
luminosity. The envelope density is parameterized in terms of the envelope accretion rate. 
We further explored those parameters using a fixed distance of 120 pc to find the best-fit 
model based on the initial model parameters obtained from the online tool (Robitaille’s 
model ID : 3005892). In our modeling, fixed model parameters among those derived from 
Robitaille’s online fitting tool are listed in Table |5j 


Robitaille et ah 


We construct the density structures of _tlny three 


envelope, and outflow cavity, following 


Robitaille et al. 


components consisting of disk, 


(2006). 


The density structure for the envelope can be described as 


p(r,6) = 




■( 


r Y 3/2 ( 1 ■ M- 1/2 ^ , 2/^BcV 1 
/ V /in / v Un r ) 


( 2 ) 


47 r(GM*i?3)i/2 \R C / v p 0 s \p 0 
where M env is the envelope accretion rate, M* is the stellar mass, and R c is the centrifugal 
radius which is usually associated with an outer radius of the disk. This centrifugal radius 


determines the flattening of envelope structu re, u is cos 6, and un is cosine o 


a streamline of infalling particles as r —y oo ([Ulr ich 
solved by the equation of parabolic motion: 


1976 


Terebev et ah 


the angle of 


1984). It can be 


Ro + — 1 ) _ ~ 


R r 


(3) 


3 See http://www.astro.wisc.edu/protostars 
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by 


The disk density p r ofile i s assumed to be flared w i th an accret ion motion characterize d 


Shakura fc Sunvaevl f 19731) . iLvnden-Bell fc Pringlel ( 1974 ). and 


Hartmann et ahl (J,1998), 


p{w,z) = Poll - 



R*\ a 


VO 


exp < — 


2 Vh\ 


(4) 


where R* is the stellar radius, and p 0 is the normalization constant calculated from 
the integral of the density over the whole disk that corresponds to Mai S k- £7 and z are the 
radius and height in the cylindrical coordinates. The disk scale height, h, is defined as h oc 
vo 4, where /3 is the disk flaring factor. The radial density exponent, a, is obtained from a 

— ft + 1 . 

We also include the outflow cavity in the model described by 

z C av = cm 1 ' 5 . (5) 


c is max/(-R c ,max tan 0 cav ) L5 , where 7? e max is the outer radius of the envelope, and # cav is 

the cavity opening angle. 


We compared the model to the photometry in Figure 0 We also obtained the binned 
IRS and PACS fluxes measured at 9, 12, 18, 25, 70, 120, and 160 pm from IRS and PACS 
spectra for more reliable fits. The best-fit is determined by the following logarithmic 
deviation, where 


N 


R = jy>iMFx,JFx,,m)\]/N. 


( 6 ) 


1=1 


F \ u o is the observed fluxes at each wavelength, F \ um is the model fluxes at each wavelength, 


N is t 
point 


le number of data points, and 00 i is the inverse of the fractional uncertainty in each 


Fischer et ah 


2012 ). 


In Figure dU we present the best-fit model for the deconvolved SED of IRS1. It has R 
3.1 (models with R < 4 provide acceptable fits within observational uncertainty). We 
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also projected the sub-mm fluxes at observed aperture sizes. The model overestimate of the 
flux compared with the data from SHARC-II 350 /mi is acc epta ble be cause the observation 


was insensitive to extended emission such as GSS30-IRS1 (Wu et al 


2007 b The model 


has a total luminosity of L tot = 10 ± 1 L & , an envelope mass of M env = 0.23 M e , and an 
envelope accretion rate of M env = 3 x 10” 6 M 0 yr -1 . The total luminosity derived from 
the best-fit model is greater than the observed Tboi of 8.5 L & , indicating that sub-mm data 


does not trace the total extent of GSS30- 


derived from sub-mm fluxes (Enoch et al 


RSI. The envelope mass falls within the mass 


2009). The envelope accretion rate derived from 


the best-fit model is the same as that found in Robitaille’s online fitting tool. The cavity 
opening angle is 5.2 ° and the inclination angle is 41 °, where 0 ° is face-on. The interstellar 
extinction a t V ba nd is 12.19, which is consistent with the mean extinctions derived fro m 


Evans et al. 


(120091) . When compared to the ID envelope model (iKristensen et al.l 2012), 


this SED model much better describes the NIR fluxes from the disks. 


After modeling the deconvolved IRS1 spectrum, we now model the PACS total emission 
over 5x5 spaxels after subtracting deconvolved emission from IRS2 and IRS3 (hereafter 
PACS total fluxes) in order to study whether the PACS total fluxes can be explained solely 
by the internal heating by adjusting the model parameters of the best-fit model of IRS1. 
In order to fit the PACS total fluxes (Figure [5l blue lines), models need to produce higher 
fluxes at FIR wavelengths and a smaller slope at longer wavelengths compared to the model 
fitting of the deconvolved SED of IRS1. When we increase the total internal luminosity in 
the model, the flux in FIR rises, but the fluxes at NIR are also enhanced (Figure USE). If 
the outer radius of the envelope increases, the slope of SED at long wavelength decreases 
(Figure fl9b). With higher envelope mass infall rates that provide high opacity, the NIR 
flux decreases because the envelope blocks the photons from the disk (Figure fl9b). (Here, 
the parameter of “mass infall rate” in the model really controls the opacity of the envelope.) 
The opening angle of the outflow cavity affects the SED only at a wavelength shorter than 
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about 100 /irn; fluxes at A < 3 /iin increase with opening angle, while fluxes at 3 /ini < A < 
100 /irn decrease with opening angle (Figure flUll) . 

None of these models are a good fit to the total FIR derived from PACS, suggesting 
that other input parameters are needed. We also tried to fit the full SED using Robitaillc’s 
online SED fitting tool in order to confirm the above results because Robitaillc’s model 
deals with only the internal heating. The minimum chi-square model fails to fit the observed 
fluxes beyond 100 /irn wavelengths. We note that before the outflow cavity is carved out, 
the envelope density structure of our model is assumed to be spherically symmetric in 
spite of two other YSOs located closely to IRS1; three YSOs are possibly gravitationally 
interacting to deviate the density structure from spherical symmetry. Geometrical effects, 
such as aspherical distributions may play a role, but the 70 /irn emission from internal 
heating always arises on small scales. Only external heating can provide the extended 
excess emission that we observe. 


External Heating 


The possibility of an enhanced radiation held in the Oph cloud compared to the 


standard rac 
by 


iation held in the solar neighbo rhood has been previously studied in depth 


Yui et al. 

(1993 

) and 

Liseau et al. 

(1999) 


(119991) . They suggested that the Oph region might 


be filled with strong UV radiation on a large scale, indicative of diffuse PDRs over the 
entire Oph region. 


Yui et ah 


(119931) found that the UV radiation extends a few parsecs 


from the nearby B type star over the whole cloud. 


Lindberg & Jprgensen (2012) reported 


that low-mass protostars cannot heat the surrounding envelope up to 30 K, but an external 
radiation held can explain the high temperature of 40 K that we found in the residual 
spectrum. In addition, our deconvolution results show that the huxes measured over the 
PACS 25 spaxels greatly exceed the sum of each deconvolved hux of three YSOs, suggesting 
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that external heating may contribute to the total flux density of GSS30-IRS1. Therefore, 
we discuss how the external radiation field has an effect on the PACS total IRS1 fluxes. 

The modeling of the SED with both internal and external heating was done with the 
same physical parameters obtained from the model with only internal heating. We also 
modeled the PACS fluxes in their observed apertures, because the PACS FOV cannot trace 
the outer region of the source where the external heating is more dominant than internal 
heating. 

First, we included an external radiation field of 20 Co, expected to be the minimum 
Co from [C 11] emission (Sect. 14.21) . The best-fit SED model with 20 Co is presented in 
Figure [20] (blue line). However, the model fails to fit the observed fluxes in the PACS 
apertures (blue squares). Furthermore, we adopted the simple blackbody spectrum of a B2 
V star as an external radiation field, and we set its position using the projected distance 
between the B2 V star and GSS30-IRS1, approximately 0.48 pc given a distance of 120 pc 


to the Oph c. 


(Liseau et al. 


oud. If the effective temperature of the B2 V star is taken to be ~ 22,000 K 


1999|), the star increases the radiation field by a factor of 226 Cq at a distance 


of 0.48 pc. In this model, the total fluxes of the best-fit model correspond to the observed 
fluxes, but the model still cannot reproduce the fluxes in the PACS apertures (orange 
squares in Figure [201) . According to our models, the UV flux must be enhanced more than 
1000 Co in order to fit the PACS data (red squares in Figure [20]). If the GSS30-IRS1 and 
the B2 V star are not in the same plane of sky, the 226 Co, which is calculated at the 
projected distance of 0.48 pc, would give us the maximum C 0 . Therefore, the value of 1000 
C 0 is unrealistic. This test suggests that if heating only by the high energy photons is 
considered, the B2 V star cannot be the heating source. 


Finally, we adopted the spectrum of ISRF from Ev ans et a l 


external heating effect. They combined the spectra from 


(2001) to explore the 


Black ( 1994 ) for A > 0.36 /iin 
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and Draine (11978!) for A < 0.36 /irn (Figure [2TJ the Black-Draine field). After modeling 
the SED with the ISRF enhanced by factors of 5 to 500, we find that the model with the 
1SRF enhanced by a factor of 130 can reproduce the observed PACS total fluxes well. If 
we calculate the Go value at the FUV range using this enhanced ISRF, it is 201, which is 
much smaller than the required enhancement (1000 Go) when the spectrum of the B2 V 
star is adopted. Nevertheless, the ISRF enhanced by a factor of 130 can fit the observations 
because the ISRF spectrum has much higher intensities at wavelengths greater than the 
wavelength of intensity peak of the B2 V star as seen in Figure 1221 

The best-fit SED model with the ISRF enhanced by a factor of 130 is presented in 
Figure [231 Our model photometry at 70 /irn fits the observed total fluxes well, but there 
is a disagreement at 160 /irn. It is possible that this flux deficit at 160 /irn is caused by 
a larger structure surrounding GSS30-IRS1. In the SCUBA dust continuum map at 850 
/irn (Figure HJi), the emission is extended beyond the PACS FOV and its peak is closer to 
IRS3 than IRS1. In addition, the 1.3 mm continuum flux dens ity of IRS3 after background 


subtraction is greater than IRS1 by a factor of 2 ( IMotte et ah 


1998). In order to see how 


the core structure can influence the fluxes at the PACS longer wavelength, we increased the 
outer radius of the envelope by a factor of 2 compared to the best-fit model with a fixed 
constant density. Our simple test confirmed that cold dust material in the larger structure 
leads to the increase of flux at 160 //rri, although the NIR flux decreased due to the added 
extinction. 


The dust temperature profile derived from three different models is presented in 
Figure [24]. In the model without external heating, the temperature drops below 20 K 
beyond 4000 AU, indicating that internal heating cannot produce a temperature of 40 
K, which we found from the residual spectrum, at large radii. However, the same model, 
with the interstellar radiation field enhanced by a factor of 130, can produce such high 
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temperature at the outer envelope radius. Without a central heating source, the external 
heating model cannot explain the high temperature at the inner part of envelope (< 2000 
AU). 


Our various models show that dust continuum emission from the PACS 25 spaxels 
cannot be produced solely by the internal heating from the central object or UV fluxes from 
the B2 V star, but it can be explained by the enhanced Black-Draine interstellar radiation 
held, which has intensities increasing with wavelength, indicative of the importance of low 
energy photons for the dust heating. In order to test the contribution of low energy photons 
in our SED modeling, we enhance the interstellar radiation held only for the wavelength 
range from 0.01 /ini to a certain boundary wavelength; except for this wavelength range, no 
radiation is assumed to exist. In Figure 1251 each boundary of upper wavelength is presented 
with different colors. Our test supports that low energy photons (0.25 eV < hv < 6 eV; 
0.2 /irn < A < 5 /irn) as well as the UV photons (6 eV < hv < 13.6 eV; 0.1 /irn < A < 
0.2 /irn) play a role in heating dust material in the envelope. (However, the photons with 
energies lower than 0.25 eV (A > 5 /irn) do not affect the dust heating significantly as seen 
with the red dashed line in Figure 1251 ) Therefore, stars with spectral types later than B 
can contribute to the external heating of YSOs. On the other hand, the [C II] line emission 
is predominantly produced by UV photons from the B2 V star since carbon is ionized by 
UV photons, specifically. 


Green et ah 


(120131 ) presented the ratios of flux densities extracted from the PACS 3 


x 3 spaxels and the PACS central spaxel for the 15 well-centered sources including Eve 
sources (GSS30, VLA1623, Elias29, WL12, and IRS63) in the Oph cloud. According to 
their results, the ratios in the Oph sources (except for WL12 and IRS63) greatly exceed 
the PACS/PSF correction. Therefore, the continuum models of sources in the Oph cloud 
(or other externally heated sources) need more careful modeling to avoid misfitting the 





observed data, especially when their spatial resolutions are poor. 


6. SUMMARY 

We analyzed a Class I YSO, GSS30-IRS1, observed with the Herschelj PACS, as a part 
of the DIGIT key program. We investigated its physical conditions based on the PACS 
observations by analysis of molecular and atomic lines and modeling the dust continuum 
emission using the 3D dust continuum radiative transfer code (RADMC-3D). 

1. The three YSOs around GSS30 (IRS1, IRS2, and IRS3) are decomposed from the 
total FIR continuum flux across the full PACS array (5x5 spaxels). The deconvolved 
fluxes of IRS1 are consistent with the PACS central flux, corrected for the PACS/PSF. At 
longer PACS wavelengths, the deconvolved IRS1 flux is comparable to the deconvolved 
IRS3 flux. 

2 . Almost 75 molecular and atomic lines are detected: the ladders of CO from J = 
14-13 to 40-39, 38 H 2 0 lines from 100 K to 1500 K of excitation energy, 14 transitions of 
OH lines in two different ladders, and two atomic [O I] lines at 63 and 145 /irn. The [C II] 
158 /.irn emission line, which is known as a tracer of the PDRs, is also detected. 

3. All detected emission lines peak in the central spaxel and are spatially compact, 
except for [O I] and [C II]. The [O I] emission is extended along the NE-SW direction, which 
is similar to the outflow direction reported in previous studies. However, the [C II] emission 
shows no correlation with any other line emission. One special feature of GSS30-IRS1 is 
that the continuum is extended toward the northeast direction, beyond the PACS/PSF, 
indicative of anisotropic external heating. 

4. The rotational diagrams and non-LTE LVG models indicate similar molecular gas 
properties to other YSOs observed with Herschel, indicating that any external radiation 
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field in GSS30 does not have a great effect on the excitation conditions. 

5. The [C II] emission provides information on the FUV radiation field around GSS30 
(3 ~ 20 Go)- Although the [C II] line can be excited by the UV photons from the nearby B 
type star, dust continuum emission over all PACS spaxels cannot be explained solely by a 
FUV radiation field. 

6. Our best-fit model of the PACS SED supports the idea that the enhanced interstellar 
radiation held is significant in GSS30. In this model, the dust temperature at the outer 
envelope radius corresponds to the temperature found in the residual spectrum after the 
deconvolution method. Thus, both internal and external heating must be considered to 
better constrain the physical structures of the envelope. Advanced dust radiative transfer 
models are required to understand the complex envelope morphology around GSS30, as 
well as those around other sources exposed to significant external radiation fields. 
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Table 1. Detected Line Fluxes of GSS30-IRS1 


Species Transition 

E U p 

A 

Flux a 


(K) 

H 

(10“ 18 W m" 2 ) 


40-39 

4512.67 

65.69 

130 ± 27 

37-36 

3871.69 

70.91 

73 ± 10 

36-35 

3668.78 

72.84 

131 ± 21 

35-34 

3471.27 

74.89 

92 ± 9 

34-33 

3279.15 

77.06 

146 ± 13 

33-32 

3092.45 

79.36 

152 ± 17 

32-31 

2911.15 

81.81 

155 ± 21 

31-30 

2735.28 

84.41 

348 ± 14 

30-29 

2564.83 

87.19 

194 ± 19 

29-28 

2399.82 

90.16 

188 ± 15 

28-27 

2240.24 

93.35 

230 ± 13 

25-24 

1794.23 

104.44 

196 ± 8 

24-23 

1656.47 

108.76 

232 ± 14 

23-22 

1524.19 

113.46 

661 ± 24 

22-21 

1397.38 

118.58 

282 ± 8 

21-20 

1276.05 

124.19 

306 ± 7 

20-19 

1160.20 

130.37 

285 ± 6 

19-18 

1049.84 

137.20 

279 ± 6 

18-17 

944.970 

144.78 

290 ± 13 

17-16 

845.590 

153.27 

317 ± 10 

16-15 

751.720 

162.81 

327 ± 11 

15-14 

663.350 

173.63 

320 ± 11 

14-13 

580.490 

186.00 

416 ± 13 
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Table 1—Continued 


Species 

Transition 

E up 

(K) 

A 

On) 

Flux “ 

(10“ 18 W m" 2 ) 

OH 

19 17 

2 ’ 2 2 5 2 

875.100 

55.89 

147 ± 22 



875.100 

55.95 

126 ± 40 


3 9 3 7 

2 ’ 2 2 5 2 

512.100 

65.13 

608 ± 35 



510.900 

65.28 

204 ± 30 


17 15 

2 ’ 2 2 5 2 

617.600 

71.17 

97± 6 



617.900 

71.22 

134 ± 52 


113 3 

2 ’ 2 2 5 2 

181.900 

79.12 

200 ± 6 



181.700 

79.18 

177 ± 12 


3 7 3 5 

2 ’ 2 2 ’ 2 

291.200 

84.42 

347 ± 13 



290.500 

84.60 

307 ± 20 


3 5 3 3 

2 ’ 2 2 5 2 

120.700 

119.23 

141 ± 8 



120.500 

119.44 

190 ± 30 


13 11 

2 ’ 2 2 ’ 2 

270.200 

163.12 

58 ± 9 



269.800 

163.40 

67 ± 10 

p-H 2 0 

431-322 

552.300 

56.33 

299 ± 26 


4 2 2“3i3 

454.300 

57.64 

309 ± 42 


7 2 6~6l5 

1021.00 

59.99 

136 ± 20 


431“4o4 

552.300 

61.81 

314 ± 42 


8o8^7i7 

1070.60 

63.46 

98 ± 31 


33i~2 2 o 

410.400 

67.09 

223 ± 18 


524“4i3 

598.800 

71.07 

227 ± 17 


7l7”6o6 

843.800 

71.54 

185 ± 14 


6o6 _ 5i5 

642.700 

83.28 

248 ± 26 
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Table 1—Continued 


Species 

Transition 

E U p 

(K) 

A 

(fj.ni) 

Flux “ 

(10“ 18 W m" 2 ) 


322~211 

296.800 

89.99 

341± 20 


542-533 

877.800 

94.21 

67 ± 17 


4o4“3l3 

319.500 

125.35 

147 ± 6 


331-322 

410.400 

126.71 

32± 4 


3l3“2o2 

204.700 

138.53 

188± 6 


4l3“322 

396.400 

144.52 

61 ± 6 


431“4 2 2 

552.300 

146.92 

38 ± 7 


3 2 2“3l3 

296.800 

156.19 

137 ± 7 

o-H 2 0 

432-321 

550.400 

58.70 

452 ± 27 


b- 

o 

00 

00 

1070.70 

63.32 

249 ± 39 


7l6“625 

1013.20 

66.09 

305 ± 29 


330“221 

410.700 

66.44 

468± 16 


330-303 

410.700 

67.27 

244 ± 15 


7o7“6l6 

843.500 

71.95 

319± 13 


725-634 

1125.70 

74.95 

194 ± 28 


321“2l2 

305.300 

75.38 

589 ± 24 


423“3l2 

432.200 

78.74 

518± 34 


6l6“5o5 

643.500 

82.03 

460 ± 26 


643-634 

1088.80 

92.81 

81 ± 11 


441-432 

702.300 

94.71 

115± 13 


221 -T 10 

194.100 

108.07 

290± 16 


4 l 4 - 3 o 3 

323.500 

113.54 

661 ± 24 


432“4 2 3 

550.400 

121.72 

48 ± 7 
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Table 1—Continued 


Species 

Transition 

E U p 

(K) 

A 

w 

Flux ° 

(10~ 18 W in” 2 ) 


423“4i4 

432.200 

132.41 

93 ± 5 


5i4”5 05 

574.700 

134.93 

54± 8 


33o~3 2 i 

410.700 

136.50 

100 ± 9 


3o3 _ 2i 2 

196.800 

174.63 

155 ± 6 


2 i2 _ loi 

114.400 

179.53 

179 ± 6 


2 2 i-2 12 

194.100 

180.49 

68 ± 6 

[0 I] 

3 Pi- 3 P 2 

227.712 

63.18 

6842 b ± 165 

[0 I] 

3 Po- 3 Pi 

326.579 

145.53 

382 fc ± 21 

[C II] 

2 P3/2 -2 P1/2 

91.210 

157.74 

586 c ± 15 


Note. - The CO & OH (84 /im) and CO & o-H 2 0 (113 /mi) 
emission lines are excluded due to line blending. All lines within 
the wavelengths from 95 /im to 102 /im are not included and are 
unreliable because of high noise levels. 

a All line fluxes are extracted from the central spaxel, then 
they are corrected by the PACS/PSF. 

b Flux measured from 15 spaxels of one nod position. 

c Flux measured from 25 spaxels of one nod position. 
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Table 2. Summary of the Rotational Diagram Results 


Species 

Component 

Trot 

(K) 

■ARmolecule) “ 
(10 46 ) 

CO 6 

warm 

359 ± 4 

373 ± 12 


hot 

783 ± 32 

77 ± 12 

h 2 o 

Para 

171 ± 5 

0.17 ± 0.006 


Ortho 

197 ± 2 

0.24 ± 0.006 

OH 

2 n 3 /2- 2 n 3 /2 

193 ± 7 

0.13 ± 0.008 


2 ni/ 2 - 2 n 1/2 

128 ± 5 

0.60 ± 0.10 


a The total number of molecules. 


b CO ladders are fitted by two rotational temperatures 
with a break-point at E u p = 1800 K. 
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Table 3. Summary of the best-fit LVG model parameters 




One Component 


Power Law 


Species 

Tkin 

n(H 2 ) 

N (molecule) / Av 

x 2 

power-law index, b 

n(H 2 ) 

x 2 


(K) 

(cm -3 ) 

(cm -2 km _1 s) 



(cm -3 ) 


CO 

5000 

3.2 X 10 4 

10 15 

18.3 

2.6 

1.8 x 10 6 

8.4 

h 2 o 

500 

10 7 

3.2 x 10 16 

40.4 

2.2 

10 13 

189.1 

OH 

175“ 

1.3 x 10 9 

1.8 x 10 15 

8.1 





125 6 

1.0 x 10 8 

1.0 x 10 15 

15.5 





a The best-fit model without IR-pumping effect. 
b The best-fit model with IR-pumping effect. 
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Tablc 4. Line luminosities in FIR a 


T b 

-^Continuum 

T c T d 

-^FIR -^mol 

L [o i] 

Leo 

^h 2 o 

^OH 

530 

1.00 0.68 

0.32 

0.21 

0.36 

0.11 


a In units of 1CT 2 L©. 

b Total FIR luminosity from the whole PACS spaxels. 

c Total FIR line luminosity, Tfir = Tco + -^h 2 o + -^oh + 
L[oi]- 

d Molecular line luminosity, T mo i = Tco + Th 2 o + Tqh- 







Tabic 5. 


Model parameters for the disk and envelope model 


Parameter 

Description 

Unit 

Value 

M* 

Stellar mass 

Mq 

0.56 

R* 

Stellar radius 

Rq 

8.05 

M env 

Envelope mass 

Mq 

0.23 

Me nv 

Envelope accretion rate 

10 -6 Mq yr _1 

3.00 “ 

•^e,max 

Envelope outer radius 

AU 

8472 

Mdisk 

Disk mass 

1 

o 

9.12 

■^d,max 

Disk outer radius (Centrifugal radius) 

AU 

17.40 

hioo 

Disk scale height at 100 AU 

AU 

6.11 

a 

Disk radial density exponent 


2.08 

p 

Disk scale height exponent 


1.08 

Pcav 

Cavity opening angle 

O 

5.21 

A v 

Extinction 


12.19 

Ltot 

Total luminosity 

L© 

10.00 ° 

^incl 

Inclination 

o 

41.41 


Note. — All parameters, except for the envelope accretion rate and the total 
luminosity, are obtained from Robitaille’s online fitting tool (model ID : 3005892). 

a The best-fit model parameters derived from our explored modeling using a fixed 


distance of 120 pc. 
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Fig. 1.— Spitzer -IRAC 2 (4.5 /irn) image (a) and Herschel /PACS (70 pm) image (b) of 
GSS30. The purple and red crosses represent the positions of PACS pixels. The spaxels 
where we used for measuring the [O 1] fluxes are represented as red crosses. The green 
circles show the locations of YSOs. The orange contours for the SCUBA 850 pm data are 
overplotted in the right panel. Contour levels are increasing in 10, 20, 30, 40, 50, 60, 70, 80, 
and 90 % of the peak flux. 
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20 10 0 -10 -20 20 10 0 -10 -20 20 10 0 -10 -20 

RA offset (arcsec) 


Fig. 2.— Contour maps of selected line flux (red solid line) and local continuum (black 
solid line). The positions of PACS pixels are represented as the black crosses. In each map, 
contour levels are 30, 60, and 90 % of the peak flux of each line. The gray filled circles 
indicate the beam at each wavelength. The wavelength of each line transition is presented 
in the bottom of each box. The local continuum bands are 59-61, 63-65, 88-90, 177-179, 
145-147, and 115-117 /ini from the top left in the clockwise direction. The extended line 
structure seen in the SW of GSS30-IRS1 is caused by the outflow of VLA1623 (Sect. 13. 3p . 
Unlike other molecular lines, the [O I] line is extended toward northeast direction. All local 
continuum are more extended than lines, except for [O I]. 
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A (/xm) 


Fig. 3.— The spectra of the deconvolved IRS1 (red solid line), IRS2 (purple solid line), and 
IRS3 (blue solid line) YSOs in GSS30. Fluxes are measured from the central spaxel and 
corrected for PACS/PSF, represented as black solid lines. 
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Fig. 4.— The SED of the continuum residual in PACS wavelength range (black solid line) 
and fitting of black-body dust temperature of 40 K (red solid line). 
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Fig. 5.— SED with Herschel/ PACS, 2MASS (JHK; filled circles), Spitzer -IRAC (3.6, 4.5, 
5.8, 80 /im; fille d triangles), and S pitzer -IRS (10-36 /xm; black solid li ne) from the “c2d ” 


legacy program (jLahuis et al. 


2006) and the “IRS_Disk” GTO program ( Furlan et al. 2006). 


The red and blue PACS spectra have been extracted differently; the red spectrum represents 
the deconvolved IRS1 SED while the blue spectrum shows the SED over the whole 25 PACS 
spaxels after subtracting the deconvolved fluxes from IRS2 and IRS3. The Sub-nun data are 


van Kem pen et al. 


2009 ). and IRAM (1.3 mm in 11" beam; 


9" beam, 85( 

/irn 

Motte et al. 

1998) 
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Fig. 6.— Detected CO lines from GSS30-IRS1. Lines are extracted from the central spaxel. 
The wavelength of each transition is marked: CO (red), p-H 2 0 (purple), o-H 2 0 (blue), and 
OH (green). The rotational transition level of each line is listed on top of each box. 
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Fig. 7.— The same as Figure [6] but for p-H 2 0 lines. 
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Fig. 8.— The same as Figure [6] and [7] but for 0 -H 2 O lines. 
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.— The same as Figure [6] 0 and [8] but for OH lines. 
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Fig. 10.— The whole PACS spectrum map of [O I] 63 /jrn line. The [O I] lines southwest 
from the central spaxel are affected by the outflow of VLA1623. 
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Fig. 11.— PACS spectral maps of [0 I] 63 /in i, where each nod position spectrum is sub¬ 
tracted. 
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Fig. 12.— The same as Figure fill but for [O I] 145 /ini line. 
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Fig. 13.— The same as Figure |TT] and [T2] but for [C II] 158 //m line. The [C II] emission 
line is detected at only one nod observation. The other nod position was contaminated by 
the strong [C II] emission, resulting in the deep absorption feature in the source spectrum 
as seen in (a). 
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Fig. 14.— Contour map of [C II] 158 pm line flux (red solid line) on top of the local 
continuum emission at 151 pm (gray scale). The line flux was calculated from Figure fl3b. 
The positions of PACS pixels are indicated as the black crosses. The gray filled circle at 
the right bottom shows the PACS beam size at 158 pm. Contour levels and gray scales are 
increasing in 30, 60, and 90 % of the peak flux. 
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Fig. 15.— The line ratio of [0 I] 63 /irn and 145 /irn as a function of shock condition. 
Triangles are for n(H 2 ) = 2 x 10 4 cm -3 and squares are for n(H 2 ) = 2 x 10 5 cm ~ 3 . Black 
solid and dashed lines indicate C- and J-shock models, respectively. Red colors represent 
the flux ratio in GSS30-IRS1 and the vertical lines represent the uncertainties. 
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[0 I] 63 /i,m [0 I] 145 /xm 




Fig. 16.— The same as Figure [15] but for the average intensity of [0 I] 63 /mi and 145 /jrn 
lines. 
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Fig. 17.— The total line luminosity ratios toward the embedded sources. Left: the ratio 
between Loh and Lh 2 o- Right: the ratio of L [o i] and Th 2 o- GSS30-IRS1 is presented as the 
red crosses. 
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Fig. 18.— The best-fit SED model (red solid line) for the Herschel /PACS deconvolved 
spectrum of IRS1. Data symbols are the same as in Figure EJ Blue squares represent the 
model fluxes with the actual aperture sizes. 
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Fig. 19.— Exploration of parameter space: (a) total luminosities, (b) outer radius of enve¬ 
lope, (c) envelope mass infall rate, and (d) outflow cavity angle. The different line colors in 
each box represent the SEDs calculated with different values of a given parameter. No ad¬ 
justment of these parameters can fit the total PACS fluxes after subtracting the deconvolved 
spectra of IRS2 and IRS3. 
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Fig. 20.— Three different best-fit models (blue, orange, and red solid lines) for the total 



fluxes after subtracting the deconvolved spectra of IRS2 and IRS3. The symbols for the 
observational data are the same as in Figure [5] except for the Herschelj PACS. Squares 
represent the fluxes modeled with actual aperture sizes. 
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Fig. 21.— Plot of the standarc 


interstellar radiation field used in the dust continuum mod¬ 


eling. It is the combination of Black ( 19941 ) (for A > 0.36 /jrn) and Draine ( 1978 ) (for A < 
0.36 /jrn.) 
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Fig. 22.— The spectrum of the interstellar radiation field enhanced by a factor of 130 (black). 
For comparison, we also present the radiation field of the B2 V star at the distance of 0.48 
pc from the GSS30-IRS1 (red). 
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Fig. 23.— The best-fit SED model (red solid line) with the interstellar radiation held en¬ 
hanced by a factor of 130 for the total Herschel /PACS fluxes after subtracting the decon¬ 
volved fluxes by IRS2 and IRS3. The symbols for the observational data are the same as in 
Figure [5] except for the Herschel/ PACS. Blue squares represent the model fluxes with the 
actual aperture sizes. 
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Fig. 24.— The dust temperature profile as a function of radius. The black solid line is the 
model including only he internal heating by IRS1. The same model but with the external 
heating by the enhanced interstellar radiation field is indicated as a red solid line. The black 
dashed line is represents the model without the central source but including the external 
radiation field. 
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Fig. 25.— The SED models (colored solid and dashed lines) with the external heating by 
the ISRF for the total Herschel /PACS fluxes after subtracting the deconvolved fluxes by 
IRS2 and IRS3 (black solid line). In order to test the contribution of low energy photons, 
we set the wavelength of the ISRF differently; the interstellar radiation held enhanced by a 
factor of 130 is used in the wavelength range from 0.01 /iin to a certain wavelength, and the 
ISRF is set to be 0 in the rest part of wavelengths. Different colors represent different upper 
boundaries of wavelength of the ISRF, as marked in the upper right of the box. Squares 
represent the fluxes at 70 /iru modeled with actual aperture sizes. 











